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Abstract: We show that the B — * X s ^y photon energy (Ey) spectrum can be reli- 
ably computed by resummed perturbation theory. Our calculation is based on Dressed 
Gluon Exponentiation (DGE) incorporating Sudakov and renormalon resummation. It is 
shown that the resummed spectrum does not have the perturbative support properties: it 
smoothly extends to the non-perturbative region Ey > m/2, where m is the quark pole 
mass, and tends to zero near the physical endpoint. The calculation of the Sudakov factor, 
which determines the shape of the spectrum in the peak region, as well as that of the 
pole mass, which sets the energy scale, are performed using Principal- Value Borel summa- 
tion. By using the same prescription in both, the cancellation of the leading renormalon 
ambiguity is respected. Furthermore, in computing the Sudakov exponent we go beyond 
the formal next -to-next-to-leading logarithmic accuracy using the large-order asymptotic 
behavior of the series, which is accurately determined from the relation with the pole mass. 
Upon matching the resummed result with the next-to-leading order expression we com- 
pute the spectrum, obtain its moments as a function of a minimum photon energy cut, 
analyze sources of uncertainty and show that our predictions are in good agreement with 
Belle data. 
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1. Introduction 

Inclusive measurements of radiative, B — ► X s j, and semileptonic decays in the B factories 
have a great potential in putting stringent constraints on short-distance physics beyond 
the Standard Model and in providing an accurate determination of the CKM parameters. 
However, exploiting this potential crucially depends on our ability to extrapolate from the 
experimentally accessible kinematic domain to the full phase space. This requires precise 
theoretical predictions for inclusive decay spectra. 

The main obstacle in the QCD calculation of inclusive decay spectra [1-15] - as 
opposed to total rates — is their sensitivity to the momentum distribution of the heavy 
quark in the meson [1,2], which is defined by 

/(*,M; M ) = J°° ^-e- izP+ y~ {B(P)\ 9(y)* v (0, y) 7+ *(0) \B{P)) ^ ; P 2 = M 2 , 

(1.1) 
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where M is the meson mass, \x is the renormalization scale of the operator, y is a lightlike 
vector in the "— " direction and & y (0,y) is a path-ordered exponential in this direction 
connecting the points y and 0. Being a property of the bound state, this distribution 
is of course non-perturbative. Although it can be analyzed using the Operator Product 
Expansion (OPE) and written as an infinite sum over forward matrix elements of local 
operators, quantitative information on this distribution is very limited. 

The quark distribution function is directly relevant in the experimentally most accessi- 
ble region near the endpoint of the spectrum, where the invariant mass of the hadronic jet 
is small. In B — ► X s j decays this corresponds to the photon energy E 1 (measured in the 
B rest frame) being close to its maximal value M/2. This limit is inherently important as 
the leading-order partonic process corresponds to a photon and an s-quark recoiling back 
to back, with Ey = to/2 (to is the quark pole mass). This results in a 6(1 — x) distribution, 
where x = 2E 1 /m. This distribution is smeared by perturbative and non-perturbative 
effects but it still peaks near x = 1. Because of the singular nature of this distribution, 
it is convenient to consider the photon-energy moments. The perturbative expansion of 
these moments is well defined to any order in perturbation theory, but it is dominated 
by Sudakov logs, lniV (N is the moment index) and it therefore requires resummation. 
The dominant non-perturbative contributions, which are associated with the momentum 
distribution function of the heavy quark in the meson, appear as powers of NA/m. This 
parametric enhancement of perturbative and non-perturbative contributions at large N 
is important even if one considers only the first few moments that are measured experi- 
mentally: the information encoded in high moments is absolutely essential to recover the 
correct dependence of the partial decay rate on the minimal photon energy cut Ey > Eq. 
Such a cut is experimentally unavoidable. 

Being motivated by the OPE, most theoretical approaches to decay spectra have re- 
lied on introducing a factorization scale, which is either a hard cutoff or a dimensional- 
regularization scale, distinguishing between soft interaction that is associated with the 
momentum distribution in the meson and harder interaction that depends on the details of 
the decay process at hand. While the latter is described by perturbation theory the former 
can be parametrized given sufficient experimental constraints. This analysis is analogous 
to that of deep inelastic structure functions, but the analogy is restricted to a certain kine- 
matic domain near the endpoint; moreover, it is incomplete since decay spectra contain an 
additional source of double logarithmic corrections that is absent in structure functions. 
Unfortunately, in practice the choice of factorization scale, and even more so the proce- 
dure by which factorization (the separation between the perturbative and non-perturbative 
regimes) is implemented, strongly affects the final answer for observable quantities. A good 
example is provided by the first few moments of the photon energy in B — ► X s ^f decays 
with experimentally relevant cuts, where the factorization prescription has been the center 
of a long lasting controversy, see e.g. [14,15]. 

In our approach [16, 17] separation between perturbative and non-perturbative cor- 
rections is implemented without introducing any factorization scale. Instead, we strongly 
rely on the resummation of the perturbative expansion. The starting point is the fact 
that at the partonic level, moments of decay spectra are infrared and collinear safe; in- 
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frared sensitivity (confinement effects) shows up only through power corrections. Applying 
Dressed Gluon Exponentiation (DGE) [16, 19-22] we resum Sudakov logarithms as well as 
running-coupling effects. The conceptual step we make is to regard the perturbative cal- 
culation of the Sudakov exponent (in moment space) as an asymptotic series. This series is 
summed using Principal- Value Borel summation 1 . This amounts to parameter-free power- 
like separation between perturbative and non-perturbative corrections. Non-perturbative 
parameters controlling powers of NA/m are then defined in this prescription, and they 
loose their immediate interpretation as local matrix elements 2 . This, however, does not 
alter the fact that these corrections are associated with the quark distribution function 
nor does it undermine their universal nature. On the other hand, it typically makes them 
numerically small compared to any cutoff-based separation. 

Sudakov resummation in inclusive decay spectra [4] closely parallels threshold resum- 
mation in hard scattering processes [23-26]: perturbative corrections from real gluon emis- 
sion are singular in the soft and the collinear limits; in inclusive observables such as the 
spectral moments these singularities cancel out by virtual corrections, leaving behind finite 
but large contributions ~ C n ^ a n s ln fc N with 1 < k < 2n at any order in perturbation the- 
ory. These corrections exponentiate. Standard techniques facilitate identifying their origin 
in phase space and then systematically resumming them to a given logarithmic accuracy. 
In inclusive decays Sudakov logarithms are associated with two independent subprocesses 
(see Refs. [4, 16], and Fig. 1 in Ref. [17]); each of them can be defined and computed to all 
orders in a process-independent manner: 

1. Sn(iti; /j,), the soft function, which is the Sudakov factor of the heavy quark distri- 
bution function, summing up radiation with momenta 0(m/N) — often referred to 
as the 'soft scale' — that influences the momentum of the heavy quark prior to its 
decay. This function is defined by considering the z — ► 1 singular terms in 

/ PT (*,m; M ) = j°° %e- izp+v ~ <&(p)l *(y)* w (0,j/) 7+ *(0) !&(?)>„ ; p 2 = m 2 , 

(1.2) 

which is the perturbative analogue of Eq. (1.1) where the external state \b(p)) is an 
on-shell heavy quark. Upon taking moments one obtains 

F^ T (m; fi) = f 1 dzf VT (z, m; ^z^ 1 = H(a s (m)) S N (m; fi) + 0(l/N), (1.3) 
J o 

where SW(m; jj) sums up the log-enhanced terms to all orders and H{a s (m)) incor- 
porates the finite terms at N — ► oo. See Ref. [30] for further details. 

2. Jjy(m,fi), the jet function, summing up radiation that is associated with an unre- 
solved final-state quark jet of invariant mass squared 0{m? /N). This function is 
directly related to the large- x limit of deep inelastic structure functions. See Ref. [32] 
for further details. 

1 \t was shown in Ref. [18] that Principal- Value Borel summation is equivalent in principle to a hard 
cutoff on some Euclidean momentum. On the other hand, it makes much less dramatic an impact on the 
distribution. In particular, it is closer to truncating the perturbative expansion at the minimal term. 

2 The locality measure in a hard cutoff based approach is the scale. Here there is no such measure. 
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The hard interaction mediating the decay is totally irrelevant for the Sudakov factor. 

Ref. [16] has generalized the large- X factorization approach beyond the perturbative 
(logarithmic) level. It has been shown that, when considered to all orders, the moments 
corresponding to each of the above subprocesses contain infrared renormalons. The quark 
distribution Sudakov factor SV has renormalon ambiguities scaling as integer powers of 
NA/m, while the jet function Jn contains ones that scale as integer powers of NA 2 /m 2 . 
This implies that certain power corrections are inherent to these subprocesses. As soon as 
these power corrections are important, Sudakov resummation ceases to be purely perturba- 
tive. Ref. [16] has shown that running-coupling effects constitute a significant contribution 
to the Sudakov exponent, and that their resummation necessarily links the calculation of 
the Sudakov factor with the parametrization of power corrections. When treated separately 
each of these ingredients is inherently ambiguous, but when combined, the ambiguities can- 
cel out. 

The leading non-perturbative corrections to inclusive decays are those associated with 
the quark distribution function, Eq. (1.1). The moments of this non-perturbative function 
can be expressed as [16]: 



where F^ T (m;[i) are the moments of the perturbative quark distribution in an on- shell 
heavy quark, Eq. (1.2), the exponential factor stands for the "binding energy" effect, and 
T sums up additional, quark-mass independent non-perturbative power corrections on the 
soft scale M/N to all orders. Power corrections on the hard scale m are neglected here. 

The exponential factor in Eq. (1.4) has an important role in computing decay spectra: 
it converts spectral moments from partonic kinematics, where moments are defined with 
respect to powers of x = 2E 1 /m, to hadronic kinematics. It is well known that the pole 
mass m has an infrared renormalon corresponding to a linear 0(A) ambiguity [27,28]. 
In Ref. [16] it was shown that this ambiguity cancels out between the exponential factor 
exp{ — (N — 1)A/M} and the Sudakov factor SV in the perturbative quark distribution 
F^ T (m;n). Thus, in the product in Eq. (1.4) one recovers an unambiguous result for the 
quark distribution in the meson. In this work we make use of this observation in two ways: 

• When computing the Sudakov factor we use the Principal- Value Borel sum. The 
same prescription is then implied for A. This is implemented here when computing 
the pole mass m from the short distance mass m^s- 

• Our approximation for the Borel function in the Sudakov factor is improved based 
on the observation that its ambiguity must cancel against that of the pole mass: the 
known large-order behavior of the relation between the pole mass and itims is used 
to fix the large-order behavior of the Sudakov factor and in this way improve its 
determination beyond what is known from fixed-order calculations. 

The main purpose of the present paper is to provide a DGE-based prediction for the 
photon-energy spectrum in B — ► X s j decays, which can be confronted with experimental 





- 4 - 



data and provide a baseline for parametrization of power corrections when the data become 
precise enough. As opposed to previous attempts to describe the spectrum, here we refrain 
from making any arbitrary parametrization of non-perturbative corrections. In Eq. (1.4) 
we use F((N — 1)A/M) = 1. We only account for the "binding-energy" effect through 
A = M — m, where the pole mass is computed from mis- Barring the uncertainty in 
modeling perturbative corrections to all orders (see below), our prediction for the spectrum 
depends only on a s and on the quark short-distance mass. The total rate depends on 
additional parameters such as Myy and the charm-quark mass, but the uncertainty in 
these parameters makes a negligible effect on the distribution. 

Recent progress in perturbative calculations [29,30] facilitates computing the Sudakov 
factor with next-to-next-to-leading logarithmic (NNLL) accuracy. We therefore begin our 
study in Sec. 2.1 by presenting the NNLL results. We investigate the convergence of the 
perturbative expansion with increasing logarithmic accuracy and confirm the prediction of 
Ref. [16] that this expansion breaks down early. Then, in Sec. 2.2 we recall the formulation 
of the Sudakov exponent as a Borel sum and review the results obtained in Ref. [16] for 
the Borel function in the large-/3o limit. In Sec. 2.3 we address one of the central issues 
of this paper, namely the construction of an approximate Borel function for the exponent 
that incorporates the exact analytic functions in the large-/3o limit on the one hand, and 
has the exact expansion coefficients in the full theory to the NNLO on the other. We 
examine in detail the sensitivity of the spectrum to the assumptions made on the Borel 
function away from the origin and incorporate the constraint on the large-order behavior 
of the exponent, which we determine from the relation with the pole mass. Sec. 3 deals 
with the matching of the resummed Sudakov factor to the full NLO perturbative result 
that incorporates contributions from all the relevant short-distance operators. In Sec. 4 we 
address non-perturbative power corrections. Here we compute the pole mass m PV from the 
short-distance mass and use it to translate the perturbative moments defined with respect 
to x = 2E 1 /m PY into the spectrum in physical photon-energy units. We end up with 
a parameter-free prediction for the spectrum, which includes all that is currently known 
about the B — ► X s ~f spectrum in perturbation theory. In Sec. 5 we present our numerical 
results for the spectrum, compute its moments as a function of the photon-energy cut and 
analyze sources of uncertainty. Finally, we compare the predictions with the available data 
from Belle. 

2. The Sudakov exponent by Dressed Gluon Exponentiation 
2.1 Sudakov resummation with fixed logarithmic accuracy 

Sudakov resummation is based on the exponentiation of logarithmically-enhanced terms 
in moment space. The spectral moments of the b — ► X s j decay occurring through the 
magnetic-operator interaction O7 (see Eq. (3.2) below) can be expressed as 




1 dr° 7 ' 



dx 




C° 7 {a s {m)) x Sud(m,A), 



(2.1) 
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where dT° r ' FT (x)/dx is the perturbative distribution in x = 2Ey/m, the superscript PT 
stands for 'perturbative', meaning in particular that the initial state is an on-shell heavy 
quark with mass m, and the bar indicates that the distribution is normalized by r^f T so 
Mn=i = 1- In Eq. (2.1) all the log-enhanced terms are resummed into the Sudakov factor 
Sud(m, TV). After the resummation has been performed, one can determine the spectrum 
from the moments by inverting the Mellin transform: 

i dr°?^(x) _ fdN _ N _ PTi07 



r° 7 - PT dx 

x total 



where the contour C extends from c — ioo to c + ioo to the right of the singularities of the 
integrand. 

The Sudakov factor Sud(m, TV) is independent of the details of the hard interaction. 
Therefore, other operators contributing to the b — ► X s j process staring at NLO share the 
same Sudakov factor. Their contribution will be taken into account in Sec. 3 and in the 
numerical analysis; for simplicity, in this section we consider only Oj. 

Based on previous analysis of large-x factorization in b — > X s j [4, 16] we know that 
the Sudakov resummation formula takes the form 3 : 



- r x Jy-l _ i r ML-x)m* j 2 . . 

Sud(m, TV) = exp <^ / dx— / ~^rM a s {u 2 ) J 

[Jo 1 ~X LV(l-x) a m 2 V V J 

+B(a s {{l-x)m 2 ))-v(a s {{l-x) 2 m 2 )^ |. (2.3) 



This all-order formula depends on three anomalous dimensions: A is the universal cusp 
anomalous dimension [53-55], which is also the large- a; limit of the quark-quark splitting 
function; B is the Sudakov anomalous dimension associated with an unresolved quark jet 
with a given invariant mass, which (in combination with A) also determines the large-x limit 
of deep inelastic structure functions [23,24,32]; and T> is the Sudakov anomalous dimension 
that is associated with the momentum distribution in an on-shell heavy quark [30]. Each 
anomalous dimension can be expanded in af s as follows: 

MS/, ,2-1 



A(a s ^)) = Y,A n (^^) , 

n=l ^ 

oo / 

n=l ^ 



afV) 



7T 

P(a s(M 2 )) = E^(^V ii ) • (2 ' 4) 

n=l ^ ' 

The first few orders in these expansions are known exactly. Higher orders are known only 
in the large-/?o limit [16]. They will be discussed in the next section. 



3 Sud(m, N) differs from Sud(m, TV) in Eq. (2.1) and in Eqs. (2.11) and (2.17) below just by finite terms 
and terms that vanish in the TV — > oo limit. These terms will eventually be discarded from the Sudakov 
factor and included in the matching coefficient Cjv(a 3 (m)). 
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The known coefficients are the following. First, the cusp anomalous dimension was 
recently computed to three-loop order [29]: 



Ax 
A 2 

A, 



= C F 



(2.5) 



36 12 ' 



5iVf 
18 



C F 



_(i cm _»)„, c , + [(£ + £ cm _ 



216 



720 



C F . 



Second, the anomalous dimension of the jet function is known to two-loop order from 
calculations in deep inelastic scattering, see e.g. Refs. [31,32]: 



3C F 



(2.6) 



3155 llvr 2 5 
W + -7^ + 2 C(3) l CA 



+ 



247 vr 2 
432 ~ 36 



N f 



C F . 



Finally, the anomalous dimension appearing in the quark distribution function in an on- 
shell heavy quark was recently computed to two-loop order [30]: 



Di = C F 
D 2 = 



(2.7) 



55 9 vr 2 
108 + 4 ^ ' ~ 12 



54 



These coefficients facilitate computing the Sudakov exponent with NNLL accuracy. 
In order to obtain a fixed-logarithmic-accuracy formula to this order we first express the 
running coupling in terms of af s (m 2 ), 



7T 



7T 



1-CT 



< s (m 2 )\ A ln(l-ff) /< s (m 2 ) 



7T 



+ 



(1-^)3 (l-a) 2 ; /3b 



/3 2 MS 



+ 



1 



where 



(1-a) 2 



/% (1 " <?? V 7T 

_ l + ln(l-o-)-ln 2 (l-o-) 



a 1 m 

A) hi— 5", 

7T ^ z 



and 



/3o 2 J 



(2.8) 
+ ••• 



(2.9) 



/? 2 MS = ^C 2 F N f + 



205 r /v + 11 /vM r 4. 2857 r3 4. 79 r n 2 1415 r 2 

1152 ^ + 576 N f j ° F + 3456 ^ + 3456^ " 3456 ° A ' 
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and then integrate the A term over // 2 . The integration over x is then done (keeping only 
logarithmically-enhanced terms to NNLL accuracy) using the general formula: 



I 



, x^ 1 - 1 ^ 

dx F 

1 — x 



f3 af s (m 2 ) ^ 1 



TT 



1-X 



TT 



dujF{uj) 



1 P , 



7E F(A)- -F'(A) [- + 1% 



TT 



2 \ /3 < S K) 



(2.10) 



+ 



TT 



where 



7T 



The resulting Sudakov factor is: 

( oo 

Sud(m,7V) = expi^ 5 „(A) 



,n=0 



7T 



n— 1 ' 



(2.11) 



where the first three coefficients <7 n (A), which sum up the logarithms to NNLL accuracy, 



4 

are : 



9oW = ^ ((1 - A) In (1 - A) - \ (1 - 2A) In (1 - 2A)) 



(2.12) 



91 (A) 



Arm/' 
A) V 

^2 

A) 2 
Ai/9i 



+- 



In (1 - A) + In (1 - 2A) ) + B 1 In (1 - A) - - L>i In (1 - 2A) 



In (1 - A) + - ln(l - 2A) 



~ In (1 - 2A) - - In (1 - 2A) 2 + - In (1 - A) 2 + In (1 - A) 



1 
2 
1 



^)=^(.2-i^2A + 



1 



-2A + 2 



+#i j E 1 



1 

+ A) 

1 



1 - A 
l : 



+ Di 7s 



2 77 
1 

1 - 2A " 



A 2 1e + 1 



1-A 1-2A 
/ln(l-A) In (1-2 A) 



1 



1 - A 



+ 



/3 2 UV A-l 2A-1 
In (1 - A) 



+ 



B 2 + 
1 



1 

A-l 



A-l 



1 1 1 

+ ! ~4 + 8A-4 ~ 2A-2 



1 ) B 1 

Az 



A-l 2A-1 
1 

2 ' 4A-2 ' 2 



2 4A-2 

Ai'yE 



D 2 



I 1 1 ln(-2A + l) N, n 

+ ^ 2A-1 J Dl 



ft 



+ 



1 / /ln(l-A) 



/3o 3 



+ 



3 lln(l-2A) 



3 



A-l 2 (A-l) 4 2 2A-1 4 (2 A-l) 



4a # 



4 An expression for (72(A) was derived a few years ago [12] although the corresponding coefficients (A3, 
D2 and B2) were not yet known. We find however that the terms proportional to A\ P2 in that expression 
(Eq. (40) there) are incorrect. 
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+ (ln(l-A)-lln(l-2A)-~ 



1 



1 



+ 



A) 4 



ln(l-2A) A 

H h 

2A-1 



2A-2 

1 In 2 (1 - A) In (1 - A) A 
2 

1 

4 



8A 
1 



.2 



A- 1 
1 



A-l 



1 In 2 (1 - 2A) 
2A - 2 + 4 2A - 1 



8A 



where the coefficients of the anomalous dimensions and the (3 function are in the MS 
scheme. They are given by Eqs. (2.5) through (2.7) and Eq. (2.9), respectively. 
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Figure 1: The moment-space Sudakov factor Sud(m, N), computed using Eq. (2.11) to LL, NLL 
and NNLL accuracy, plotted as dotdashes, dashes and dots, respectively. The curves end at N ~ 33, 
corresponding to A = ^, where the Landau singularity in g n (A) appears. 



Fig. 1 shows the Sudakov factor of Eq. (2.11) with increasing logarithmic accuracy: 
5o(A) only (LL accuracy), then go(A) and gi(X) (NLL accuracy) and finally the first three 
terms (NNLL accuracy). We observe that the convergence of the series in Eq. (2.11) is 
poor. Non-negligible differences appear already in the first few moments. These will be 
partially washed out by matching the results with the fixed-order calculation. However, 
for iV 10, where lniV starts to get large, the NNLL contribution of 52(A) is already larger 
than that of the NLL term; for N ^ 18 the NNLL contribution becomes larger than the 
leadingdog term! This clearly indicates that the perturbative expansion breaks down. It 
is known from previous studies [16, 19-22] that the asymptotic behavior of Sud(m, N) is 
controlled by infrared renormalons. In Ref. [16] it was shown, based on the resummation 
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of running-coupling effects in the large-/3o limit, that the expansion in Eq. (2.11) should 
break down early (see Fig. 2 there). Having at hand the full NNLL result we confirm that 
this is indeed so. 

The conclusions are the following: first, it becomes obvious that the logarithmic accu- 
racy criterion is, at best, insufficient. Moreover, since the poor convergence of the series 
in Eq. (2.11) is associated with running-coupling effects, these must be resummed. Fi- 
nally, since the latter brings about infrared renormalon ambiguities that scale as powers 
of (NA/m) the resolution of this problem is directly linked to the question of separation 
between perturbative and non-perturbative contributions: the calculation of the Sudakov 
exponent cannot be considered as a purely perturbative matter. Similar observations in 
different but analogous physics problems have led to the development of DGE [19-22] as 
an alternative to resummation with fixed logarithmic accuracy. 

Before moving on to DGE let us conclude our perturbative analysis by summarizing the 
state-of-the-art fixed-order and fixed-logarithmic-accuracy results. First, in Appendix A 
we derive explicit expressions for the log-enhanced terms at 0(a 2 s ). These will be useful for 
checking NNLO calculations. Finally, it is straightforward 5 to invert the Mellin transform of 
Sud(m, N) in Eq. (2.11) to obtain a resummation formula with NNLL accuracy in x space. 
The integrated spectrum with E 1 > Eq = (1 — A)m/2, 



r[6 — x sl ] 



dx 



dr°?' PT (x) 



l-A 



dx 



(2.13) 



is then: 



r[6 



*.7iS;>*(A) 



h0 7 ,pt 

total 



= C£(a a (m))- 



exp <^ g (uj) » ' ' 



i / \ i / \a„ (m) 



where 



+ I^(a;)/? 2 ^M(* 2 (1 - /^H) -*'(!- o9 ' o (u)) ) 

af(m 2 ) . . 1 

7T A 



+ RoAA), (2.14) 



Eq. (2.14) was matched to agree with the exact NLO expression upon expansion us- 
ing Eq. (A.2): 



- 1 - Or ^> (I + g - L, E + I) + 0(al) (2.15) 
incorporates finite terms in the N — > oo limit, while the additive contribution 
Ro 7 (A) = 



> 3 + U lnA+ i) A2 + U- lnALA 



Cf <M + ( a 2) (2.16) 



7T 



completes the terms that vanish at A — ► 0. Note that the currently unknown 0(a 2 s ) term 
in Eq. (2.15) influences NNLL terms in the spectrum (at 0(a^) and beyond) by mixing 



3 The general algorithm is explained in Sec. 3.4 of Ref. [19]. 
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with the leading logarithms from the expansion of the exponential function in Eq. (2.14). 
Thus, although the Sudakov factor itself is known to NNLL accuracy, complete NNLL 
accuracy of the spectrum is not yet available. 

The resulting spectrum with the Sudakov factor computed at NNLL accuracy is shown 
in Fig. 2 together with lower logarithmic accuracy results. We find that in spite of matching 
the resummed result to the NLO expression (this compensates for some of the difference 
observed in Fig. 1) differences between the three approximations to the spectrum are still 
large. Note that the three curves are not plotted beyond x ~ 0.97, where Eq. (2.14) becomes 
complex owing to the presence of the Landau singularity 6 at u> = 1/2. We observe that 
with increasing logarithmic accuracy the differential spectrum widens up and gets shifted 
to lower photon energies. While these features are intuitively expected, the instability of 
the perturbative result makes it hard to rely on this spectrum and use it as a baseline for 
parametrization of non-perturbative effects. An explicit separation criterion between the 
perturbative and non-perturbative regimes is clearly missing here. As we shall see below, 
these deficiencies are cured by DGE. 

2.2 Borel representation of the exponent and the large —/3q limit 

According to Ref. [16] the Sudakov factor in Eq. (2.1) can be expressed as the following 
Borel sum: 

Sud(m,A0 = j^jH^T^) (J^j B s (u)T(-2u) (N 2u - l) - 



Bj(u)T(-u)(N u -l) 



.(2.17) 



Here we use the scheme-invariant Borel representation [33] where T{u) is the Laplace 
transform of the 't Hooft coupling: 

with 5 = Pi/ 0q, where the first two coefficients of the (3 function are given in Eq. (2.9). 

The functions B$(u) and Bj(u) are the scheme invariant Borel representations of the 
anomalous dimensions of the soft (quark distribution) and the jet functions, respectively. 
Defining the Borel representation of the anomalous dimensions introduced in the previous 
section, namely 

A(a s (p)) = ^f r duT{u) (^) B A {u), 
Po Jo \H J 

B{a s {^)) = ^ duT(u) (^) U B B {u), 
Po Jo \M / 



6 Ref. [13] suggested to remove these singularities by resumming a set of 7r 2 terms that arise by evalu- 
ating the coupling at time-like momentum. These terms, similarly to other running-coupling effects, are 
resummed by DGE. 
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Figure 2: The differential (top) and the integrated (bottom) spectra based on the fixed- 
logarithmic-accuracy formula of Eq. (2.14), which is matched to NLO. The LL, NLL and NNLL 
accuracy results are plotted as dotdashes, dashes and dots, respectively. The three curves end at 
x ~ 0.97, where the resummed results become complex owing to the Landau singularity at lo = 1/2. 



Po Jo 



2\ u 



(2.19) 
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we have: 



Bs(u) = B A (u) - uB v (u), 
Bj{u) = B A {u) - uB B (u). 



(2.20) 



These are exact relations. 

Using the perturbative expansions in Eqs. (2.4) through (2.7), one obtains the following 
expansions for the corresponding Borel functions at small u: 



B A {u) = 1 + ( 3 4 r-2 



11 



+ 



1 



+ C 3 



2! 



+ 



(2.21) 



where c n represent contributions that are subleading in /?o; C2 and C3 are (see Appendix C 
in [30]): 

- 1 / ' 

C2 



Ca (I 
A) \3 
1 



7T 

12 



C3 = ~R 

P0 

1 

+ M 
1 



/251 7 o 

1 7T 

V288 144 



T Cs + 720^ 



c\ + 



235 11 , vr 2 

1 C3 H 

96 4 Sd 16 



32 F 



512 192 / A V 64 192 / A 128 F 



and 
B v {u) 

B B (u) 



9 + ^U C3 



3 

'4 + 



6 



72 



7T 2 




12 






( 73 


A) 


V~144 



0(u 2 



ir 2 247 / 73 5 \ C F ( 3 3 A ^ 2 

1 — I 1 A — I Fn 4 



+ 2C3 



+ 



00 



32 



<3 



(2.22) 

(2.23) 
u + 0{u 2 ). 



Upon expanding the entire square brackets in Eq. (2.17) in powers of u using Eqs. (2.20) 
through (2.23) one recovers the same log-enhanced terms presented in the previous section 
to NNLL accuracy. However, upon performing the Borel integral (barring renormalon 
singularities, see below) one can resum running-coupling effects to all orders. In the large- 
st) limit the Borel functions are given by [16,22]: 



B s {u) 
Bj(u) 



e3«(l-«) + 0(1//%), 



■ es 



1 



+ 



1-u l-u/2 



sin TTU 



TTU 



(2.24) 



Note that these are entire functions, free of singularities in the whole complex plane. In 
the full theory, i.e. for finite (3®, analytic expressions of this kind are, of course, not 
known. In Sec. 2.3 we shall construct approximations to these unknown functions based 
on their small-it expansions quoted above and some additional constraints. Because of 
the significant contribution of running-coupling effects [16, 20-22] the difference between 
performing the Borel integral (DGE) and extracting the leading logarithms to a fixed 
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logarithmic accuracy is very significant. This will be explicitly demonstrated in the next 
section. 

Writing the Sudakov factor of Eq. (2.17) in a factorized form, i.e. 



Sud(m, N) = SN(m; fx) JAr(m; fi), 

requires subtraction of logarithmic singularities: 

A 2 \" 



(2.25) 



Mm;/.)-«p{f jf x (2.26) 

B s (u)T(-2u) (JV 2 » - 1) + B A (u)\nN \, 



J/v(m; fi) = exp 



Po Jo u \m 



(2.27) 



Bj(u)T(-u) (N u -1) + 



in 



-j) B A (u)lnN 



where, as above, B^{u) is the Borel representation of the cusp anomalous dimension. This 
anomalous dimension determines the factorization scale (p) dependence of the separate 
soft and the jet functions, while Sud(m, N) is independent of \i: 



din S^im; /x) 



d\nJ N (m;fi) _ C F 

j; — 2 = —z-A(a a (ii))liLN. 

d m /i z po 



(2.28) 



Finally, note that the Sudakov factors Eq. (2.26) and Eq. (2.27) (and likewise Eq. (2.17)) 
have infrared renormalon singularities. As usual these singularities occur at integer (in both 
Jtv and Sn) and half integer (in Sn) values of the Borel variable u. To properly define 
these resummed Sudakov factors one needs to deform the integration contour off the real u 
axis or take the Principal- Value prescription. We shall adopt the latter, which guaranties 
that Sud(m, N) is a real-valued function, namely 



Sud(m,A0 



Sud(m,Af*) 



(2.29) 



and, in particular, that it is real for real positive N, as it is at any finite order in perturbation 
theory. By choosing this prescription we have made an explicit separation between what 
will be regarded as resummed perturbation theory and additional non-perturbative power 
corrections, which we shall discuss in Sec. 4 below. 

2.3 Renormalons beyond the large — po limit 

The all-order structure of the Sudakov exponent is summarized by Eq. (2.17). DGE makes 
use of this structure in spite of the limited knowledge of the functions B$(u) and Bj(u). 
In the large-po limit these functions, given by Eq. (2.24), are free of infrared renormalon 
singularities. Being anomalous dimensions they are expected to be free of renormalons also 
in the full theory. Similarly to other cases [19-22,32], infrared renormalons appearing in the 
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Sudakov exponent (2.17) have a very specific origin: the integration over the longitudinal 
momentum fraction z in the z — ► 1 limit. This integration gives rise to the factors T(—2u) 
and r(— u) in Eq. (2.26) and Eq. (2.27), respectively. 

Eq. (2.26) therefore implies that the soft Sudakov exponent has simple renormalon 
poles 7 at all integer and half integer values of u, except where Bs{u) vanishes. Similarly, 
according to Eq. (2.27), the jet function has simple renormalon poles at all integer values 
of u, except where Bj{u) vanishes. It turns out that in the large-/?o limit these functions 
do vanish at some of the would-be renormalon positions: Bs{u) vanishes at u = 1 while 
Bj{u) at all integers u > 3. It is not known whether the corresponding renormalons do 
appear in the full theory. 

The theoretical interest aside, the renormalon structure of the Sudakov exponent, 
and in particular, that of the soft function, is important for phenomenology. Obviously, 
the most important renormalon singularity is the one at u = \. It has been shown in 
Ref. [16] that the associated ambiguity cancels out between the soft Sudakov exponent and 
the leading non-perturbative correction A = M — m. For this cancellation to be realized 
- thus avoiding a spurious 0(NA/m) artifact — both quantities need to be computed 
as asymptotic expansions, regularizing the u = \ renormalon in the same manner. As 
explained in the introduction, in this paper we shall implement this idea. We shall use 
the Cauchy Principal- Value prescription in the calculation of A from the MS mass (see 
Sec. 4.2) and in the calculation of the Sudakov exponent. 

Power accuracy is, however, not easy to achieve. The difficulty is that in order to 
accurately compute the Principal Value of the Borel integral, say in Eq. (2.26), the renor- 
malon structure, including its overall normalization (the residue) must be known. Standard 
perturbative expansions of B$(u) as in Eq. (2.23), yield power series in u, which may not 
be reliable near u = \. In the following we will address this problem. Higher renormalon 
singularities in the soft function scaling as higher powers of NA/m are also relevant. One 
expects that these effects will mainly be important in the endpoint region. However, it is 
hard to know a priori how far from endpoint their influence extends. For the jet function 
the situation is different: the sensitivity to the functional form of Bj(u) away from the 
origin is small. This correspond to the fact that power corrections on the scale m 2 /N are 
smaller. In what follows we therefore concentrate on B$(u) and, first of all, on its value 
at u = \ which determines the corresponding renormalon residue in Eq. (2.26) and in 
Eq. (2.17). 

We first note that direct evaluation of the available NNLO expansion of Bsiu) at u = ^ 
using Eq. (2.20) with Eqs. (2.21) and (2.23) is unreliable. The apparent convergence of 
the series (in the Borel plane) at the first three orders is slow: for example with Nf = 4 
the terms of increasing powers of u are: 1, —0.0188 and —0.3415. As the expansion of an 
anomalous dimension this series is expected to converge, but since the NNLO is sizable 
while the large-order behavior is not known, it appears that such a direct evaluation of 

7 Simple poles appear only in the scheme invariant formulation of the Borel transform, where the Borel 
variable is conjugate to the logarithm of the scale. In the standard formulation, where the Borel variable 
z is conjugate to 1/A (where A = q s /3o/7t), the singularity transforms into a cut which is controlled by 
8 = /3i//3o owing to Eq. (B.13). 
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the residue cannot be accurate. 

Next, we note that the exact result of the large-/3o limit, Eq. (2.24), constrains the 
anomalous dimensions to O(u ) only, and that the corrections to B$(u) at 0{u) and 0(u 2 ), 
which are subleading in (3q, are large. Thus, a naive non-Abelianization approach, in which 
O(l/Po) terms are neglected, is not sufficiently accurate 8 . A similar conclusion was reached 
before [32] considering the case of the jet function. 

Fortunately, indirect information on B$(u) at u = \ is available owing to the exact 
cancellation [16] of the corresponding renormalon ambiguity with the one in A, or in the 
pole mass. We note that the singularity structure of the u = \ renormalon in the pole 
mass [34] — a simple pole in the scheme invariant Borel function — matches exactly the one 
of the soft Sudakov exponent, as it should. As discussed in Appendix B, the normalization 
of the pole-mass renormalon at \ is well under control. This allows us to fix the value of 
the Borel function: 

B s {u)\ u=h = f e*/ 2 . (2.30) 

The normalization constant q — see Eq. (B.6) — has been computed from the perturbative 
relation between the pole mass and the MS mass by several authors [35-38], obtaining good 
numerical convergence already at the available NNLO. In Appendix B we summarize our 
own study of this renormalon. Previous work on the subject was restricted to the standard 
Borel representation and to the MS scheme. We extend this analysis using the scheme- 
invariant formulation of the Borel transform. This provides an independent check of the 
accuracy of this calculation and facilitates the comparison with the Sudakov exponent. 

The result for the normalization of the pole-mass renormalon at ^ is shown in Fig. 13 in 
Appendix B. The close agreement between different calculational procedures based on the 
expansion of m/rnggfrnis) demonstrates the reliability of this determination. We conclude 
that the residue can be computed with ~ 2 — 3% accuracy over a wide range of (3q values, 
in agreement with Refs. [35-38]. As anticipated, the determination that relies on the soft 
Sudakov exponent is less accurate. Nevertheless, it does yield similar values. 

As explained above, in order to evaluate the Sudakov exponent of Eq. (2.17) using 
the Principal- Value prescription we must know B$(u) (and Bj{u)) as a function of u 
away from the origin. Any uncertainty would translate into uncertainty in the computed 
spectrum, and an ambiguity in the separation between perturbative and non-perturbative 
corrections. To gauge the numerical significance of this issue, let us construct a few models, 
which we generically call NNLL-DGE, that all share the same expansion in powers of u 
up to the NNLO but differ away from the origin. 

A natural possibility [19, 20, 22, 32] is to start with the analytic form of the large-/?o 
limit, Eq. (2.24), and include a multiplicative correction factor that modifies the expansion 
coefficients at 0(u) and at 0(u 2 ) by terms that are subleading in Pq, so as to match the 
exact coefficients given by Eq. (2.20) with Eq. (2.23). 

8 In general the naive non-Abelianization approach works well for perturbative expansions that are dom- 
inated by renormalons, not for anomalous dimensions. 
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At 0(u) (NLL-DGE) this was already done in Ref. [16] following Refs. [19-22] as 
follows: 



B 



NLL-DGE 



( u ) = e (f+ C2 )"(l - u) + 0(u 2 ), 



(2.31) 



and 



5 NLL-DGE (n) = I e (§+c 2 )« 



1 1 

+ 



sin TTU 



J " 2~ \l-u ' l-u/2) TTU 

Proceeding to 0{u 2 ) (NNLL-DGE) one can write: 



+ 0{u 2 ). 



(2.32) 



B< s\ u ) = es"(l - u) x exp I c 2 u + 



Bf{u) = 



es 



+ 



l-u l-u/2 



2 C A / 5 2 7 9 
C3 " C2 + ^ V 18" + 9"2 C3 
sin 7rti 



2! 



+ 0(u 3 



TTU 



(2.33) 



exp <^ c 2 u + 



2 C A /29 2 43 r \ C F / ^ 2 3 oA 



2[ 



+ 0(n 3 ) . 



It is straightforward to check that the this model coincides with the exact analytic functions 
in the large-/?o limit on the one hand, and has the exact expansion coefficients to the NNLO 
on the other. 

Let us recall that a similar exercise has already been done for the jet function in 
Ref. [32], in the context of the large-x limit of deep inelastic structure functions. Since in 
the case of decay spectra the soft function plays a dominant role, we shall adopt the model 
of Eq. (2.33) for Bj{u) and not consider here other possibilities. 

An alternative model for the soft function is given by 



Bf{u) 



et"(l 



u) 



1 + c 2 u + 



■ CU/ 5 2 7 9 
C3 + ^U 7r + 9"2 C3 



^ + 0(u*)\. (2.34) 



The main difference between B^\u) and Bg\u) is that the former inherently suppresses 
contributions from the large-ii region, which is, in any case, not well controlled. While in 
general the large-ii region is suppressed in the Borel integral (2.17) by (A 2 /m 2 ) u , at large 
N this is replaced by (N 2 A 2 /m 2 ) u . This suppression is relevant up to iV ~ m/A. Beyond 
this region the Borel integral still converges thanks to the suppression by the factor T(— 2u), 
however, the perturbative result is no more a valid approximation; the power expansion 
breaks down. 

For N ^ m/A the perturbative calculation by DGE is under control. Nevertheless, 
non-negligible contributions may still arise from the vicinity of the renormalon positions 
where the factor T(— 2u) is singular. These are power suppressed terms. Indeed, looking 
in Fig. 3 at the contributions to the Sudakov exponent from different sections along the 
Borel integration axis for N = 20, models (a) and (b) appear rather different: while model 
(a) is characterized by very small contributions beyond u ~ 1.5, in model (b) these are 
not small up to u of a few tens. This directly reflects the (exponential) suppression of the 
large-it region in B^\u) versus its enhancement in B^\u). We note, however, that the 
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Figure 3: Top: the Borel representation of the soft anomalous dimension, Bs(u), as a function 
of u for Nf = 4. Vertical lines indicate potential renormalon locations. Bottom: Contributions to 
the Principal- Value integral over u in the Sudakov exponent of Eq. (2.17) for N = 20 from the 
corresponding sections, each containing one renormalon (or none). In each figure the five curves 
corresponds to the different models described in the text: NLL-DGE as thin full line and NNLL- 
DGE: (a) as dotdashes, (6) as dashes, (c) as thick full line, and (d) as dots. 



NLL-DGE model of Eq. (2.31) does not generate significant contributions from u values 
beyond u ~ 1.5, although it has no inherent exponential suppression of the large-it region. 

To draw conclusions with regards to the net effect of the behavior of the models for 
B$(u) away from the origin we now examine the value of the Sudakov factor Sud(m, N) 
of Eq. (2.17). The result, shown in Fig. 4, is clear: models (a) and (b) as well as the 
NLL-DGE one are close up to very high moments. This means that the total effect of 
the large-u region of Bs(u), which is not under control, is moderate; in model (b) this 
is owing to cancellations between contributions from the different sections in Fig. 3. On 
the whole the DGE result is stable. Nevertheless, the differences between the models are 
non-negligible. 

Fig. 4 also shows that some difference between DGE and the conventional Sudakov- 
resummation procedure (with NNLL accuracy) develops already at low moments ./V < 10, 
and that it becomes very significant for N ~ 10 — 20. One qualitative difference is that the 
latter has a Landau singularity (N ~ 33) whereas the former does not [22]. 

A shortcoming of the models constructed so far is that they are inaccurate around 
u = i: Eq. (2.30) is not fulfilled. In order to incorporate the knowledge of Bs{u) at u = \ 
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Figure 4: The moment -space Sudakov factor Sud(m, N), as computed by DGE with perturbative 
expansions of NLL and NNLL accuracy (in the latter case the four models (a) through (d) are 
shown) and by conventional Sudakov resummation with NNLL accuracy. 



we consider the following model: 



B^'ia) - e : '."il u) x exp<jc 2 u+^ 



2 C A ( 5 2 7 9 



u 2 \ x W(u), 



(2.35) 



which, similarly to models (a) and (b), has the exact coefficients through 0{v?) and the 
correct large-/?o limit. Here, however, the additional factor 



W(u) 



is constructed such that correct value of B$(u) at u = \ would be reproduced, at least for 
the physically relevant values of Nf. To this end we set: 



wi = 1.144CU/A), W2 = -2.8CVA). 



(2.36) 



The values of Bs(u) at u = \ in the different models are shown as a function of l//3o 
in Fig. 5. Returning to Fig. 4, we observe that fixing the value at u = ^ makes a large 
effect on high moments N ^ 20. This shows that this information is relevant for the final 
spectrum. Nevertheless, the fact that differences with the other models are moderate for 
N ^ 20 is reassuring. 
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Figure 5: The value of the Borel representation of the soft anomalous dimension Bs(u) at u = h, 
which is proportional to the residue of the u — h renormalon in the pole mass, plotted as a function 
of l//3 . The exact result of the large-/3o limit is denoted by a cross. The calculation relying on 
perturbative expansion of m / '™ms( to ms) using the scheme-invariant Borel transform as determined 
using Eq. (B.20) is shown by plus signs. The four curves corresponds to the different models for 
NNLL-DGE described above: (a) dotdashes, (b) dashes, (c) full line and id) dots. Model (c) 
matches the computed value of Bs(u — |) at Nf — 4 by construction. 



Finally, let us focus on the peculiar feature of the large-/?o result for Bs{u), namely 
its vanishing at u = 1 , leading to the absence of a corresponding renormalon ambiguity. In 
the models considered so far we assumed that this property is shared by the full theory. 
This, however, has not been proven. Moreover, in Ref. [43] it has been shown that a u = 1 
renormalon appears in the kinetic-energy operator once terms that are subleading in l//3o 
are taken into account. This suggest that the same might occur in the Sudakov exponent, 
although this should be checked explicitly. In order to estimate the numerical significance 
of this issue, let us construct another model, (d), which has Bg"\u = 1) ^ 0: 

4 d) (n) =exp{Q+c 2 )u+i 

While B^"\u) does not respect the large-/?o limit result, it does have the correct NNLO 
expansion at u = 0. Returning again to Fig. 4, we observe that the effect of this modification 
is moderate, although probably not negligible. 

Finally, let us compare between the different calculations of the resummed Sudakov 
exponent described above at the level of the normalized differential spectrum. To this 



2 C A f 5 2 7 9 
C3 " C2 + ^U 7r + 9"2 Cs 



. (2.37) 



- 20 - 




Figure 6: The perturbative differential spectrum in b — > X s ^f with different approximations for 
the Sudakov factor as a function of x — 2E 1 /m 1 normalized by the total rate with E 7 > m/20. 
The different curves correspond to conventional Sudakov resummation with NNLL accuracy (thin 
dashed line) and DGE with perturbative expansions of NLL (thin full line) and NNLL accuracy 
(thick lines) where the three models (a), (c) and (d) are shown by dotdashed, full and dotted lines, 
respectively. 



end we match the resummed expression for the Sudakov factor with the full NLO result as 
explained in Sec. 3 below. Both the Borel integration in the Sudakov exponent in Eq. (2.17) 
and the inverse-Mellin integration in Eq. (2.2) are performed numerically, avoiding any 
further approximation. The results are shown in Fig. 6. Using DGE, the differences 
between different models with NLL or NNLL accuracy are moderate. On the other hand 
the result obtained by conventional NNLL Sudakov resummation is entirely different in 
spite of having the same formal logarithmic accuracy! 

It is interesting to note that the DGE result does not have the perturbative support 
properties: while the perturbative coefficients at any given order vanish for 22 7 > m/2, the 
DGE-resummed spectrum does not. It peaks for £ m/2 and then smoothly crosses the 
perturbative endpoint = m/2 and drops to zero at E~ = (m + 0(A)) /2. The resummed 
expression having different analytic properties than the coefficients is not surprising given 
that, strictly speaking, the sum does not exist. Nevertheless, it is remarkable that by 
taking the Principal- Value prescription in moment space we obtain a spectrum, which 
qualitatively corresponds to the decay of a higher-mass state. 
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3. Matching the resummed spectrum to the full NLO result 



In the previous section we discussed the QCD description of the endpoint region neglecting 
0(1/N) corrections. In order to recover the correct spectrum away from the endpoint 
region, the resummed result must be systematically matched onto the fixed-order pertur- 
bative expansion, which is available in full at NLO — see Refs. [45-47] and references 
therein. 

3.1 The NLO result 

The NLO calculation in the Standard Model is based on the effective Hamiltonian 

n eS = ^-v t *v tb J2c k o k (3.i) 

where Cu are the Wilson coefficients and are local operators of dimension 5 or 6. 
The most important operator contributing at leading order corresponds to the magnetic 
interaction, 

Or = {8LO*"b R ) F, u , (3.2) 

where ttIms is the running quark mass in the MS scheme and F^ u is the photon field strength. 
The complete basis of operators can be found for example in Refs. [45,47]. 

Sudakov logarithms originate in the universal soft and collinear limits in which the 
specific structure of the hard interaction is irrelevant. Therefore, the calculation in Ref. [16] 
and the formulae of the previous section, although derived starting with Oj, apply to all 
operators. Terms that are finite (or vanish) at large N, e.g. those of Eq. (A. 2), are different 
between the different operators. These terms will now be incorporated in the process of 
matching the Sudakov exponent to the full NLO result. 

Our treatment of the endpoint region is based on moment space. On the other hand 
contributions that are non-singular at x = 1 may be included either in moment space 
or directly in x space. Since b — > s^fg contributions associated with Os are singular at 
x — ► (when the photon gets soft) full NLO analysis in moment space is excluded. We 
therefore choose a mixed matching procedure where the dominant contribution at large 
x is taken into account in moment space, but some terms, which vanish at x — ► 1, are 
included directly in x space. 

Let us express the partonic decay rate with photon energy above some fixed energy 
cut E = (l- A)m/2 as: 



r[6 — ► x s j] Ej>Eo = r LO x 

where 9 (A > 0) is the heaviside function and 



\M v \ 2 e(A > 0) + R{A) 



(3.3) 



r LO = 



32vr 4 



V*V tb f(c^(» b ))\ (3.4) 



where m = nib, the b-quark pole mass, m^(m) is the running mass evaluated at m and, 
finally, [i b ~ m. In Eq. (3.3) M.y corresponds to the amplitude of the process b — ► 57 so 
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it begins at 0{aP s ) and includes all the purely virtual corrections at higher orders, while 
-R(A) corresponds to partonic processes with at least one gluon in the final state, and it 
therefore begins at 0(a s ) and depends on A. At NLO one has [45]: 



M 



v 




(w>) L 



ri + 7-? ° ff In — 



11(A) = 



C F a s 



7T 
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i < J 




+ 



(3.5) 



where 



(0) off 



dj — 



(w>)cf e V 



(4 0)eff (^))' 



The expressions for the LO coefficients C^ cB , the NLO coefficients C^ cS and the (com- 
plex) constants r$ as well as the anomalous dimensions 7^ off can be found e.g. in [44,45]. 
In R(A) we exhibited explicitly the terms 9 that are singular at the phase-space boundary 
A — ► (or x — ► 1). These terms must be resummed to all orders incorporating the 
cancellation of divergences between real emission and virtual corrections. All the functions 
4>ij(A) vanish at A = — see Appendix C for explicit expressions. 

3.2 The matching procedure 

Our goal is to improve the determination of T[b — ► ^s7]_E 7 >_B =(i-A)m/2 m Eq. (3.3) 
for small values of A by performing Sudakov resummation. At the same time we require 
that upon expansion the matched expression would coincide with the fixed-order result 
for any A. To perform the matching we first split the real-emission terms into two parts, 
writing 



1 



T[b — ► X sl ] Ei>Eo=xm/2 = \M V \ 



9(x< 1) + Gi(x) 



+ G 2 (x), 



(3.6) 



where both Gi^ix) begin at 0{a s ) and we require that G\{x) contains all the singular 
terms for x — ► 1 whereas G2{x) = C((l — a;) 2 ) so that it does not contribute to the 
differential rate near x = 1: 



1 



dT^jx) 
dx 



\M 



v\ 



5(1 - x) 



dGi(x) 
dx 



dG 2 (x) 
dx 



(3.7) 



Finally, we replace the square brackets by a matched Sudakov resummation (or DGE) 
formula in moment space while the G 2 term is left in x space, 



1 dT PT (x) 



LO 



dx 



J c— too 



PT 
N 



dG 2 (x) 
dx 



(3. 



9 These specific terms originate in O7. However, at higher orders there are similarly singular terms from 
all operators. 
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where is proportional to the Sudakov factor Sud(iV) in Eq. (2.17). Here we used the 
inverse Mellin transform formula of Eq. (2.2). 

Eq. (3.8) is, by construction, consistent with both the fixed-order expansion and the 
resummed expression in the x — ► 1 limit. From Eqs. (3.6), (3.7) and (3.8) it follows that 



M V N T = 1 + G 1 {x 



0)-/ 
Jo 



dxx 



N-l 



dGi{x) 
dx 



(3.9) 



+ 



where { } + stand for the "+" prescription (see Eq. (A. 5)). Note that the contribution to 
the N = 1 moment from the integral vanishes, so Mn=\ = 1 + G\{x = 0). We emphasize 
that, in contrast with our definition of Mn in Sec. 2, Mn=i / 1: the distribution is not 
normalized. 



3.3 Matching formulae at NLO 

This matching procedure can be applied at any order 10 . Specifically, at NLO we have: 
G 1 (x) = ^(-hn 2 (l- X )- 7 1 Hl-x)+ £ cy!fc-(A = l-s))+0(a2), 



i < j 



G 2 (x) = 



CfOs 



TV 



c ij £ ij (A = l-x)+0(a 2 s ), 



(3.10) 



i,j = 1.. 
i < j 



where 4>ij(A) = r)ij (A) + £y(A) with £ij(A) = 0(A 2 ); the explicit expressions for ijij and 
£ij appear in Appendix C. Using Eq. (3.9) with Eq. (3.10) we then get: 



= 1 + ^£. 



7T 



dxx 



N-l 



ln(l-x) 7 1 v - dr]ij(A) 

lT~i- + ^ ° ij d~A 



1 — x 



+ ^ c^(A = l)l + 0(a 2 ) 

j,j = 1..8 ^ 



i,j = 1..8 

» < i 



A=l-x 



+ 



i < j 

Performing the x-integration and extracting the single and double lnA^ terms which are 
included in the Sudakov exponent we finally obtain the matched Sudakov-resummed mo- 
ments: 



M P N T 



x Sud(m,N), 



(3.11) 



where 



1 



1 



i,j = 1. 
i < j 



7 



/(AO = " j (*(*0 + + fiW ~ 12 + 4 (lI/(Ar) + 7fi) + 2 ^ " U " 7£ J lnN 

Note that at NLO one simply has R(A = 1 — x) = Gi(x) + G2(x), however, starting at NNLO the 
relation is more complicated. 
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Figure 7: The differential spectrum according to Eq. (3.13) with M^ T computed by NNLL- 
DGE using Eq. (3.11) (where the Sudakov factor is evaluated by Eq. (2.17) using model (c) of 
Eq. (2.35)). We show the spectrum matched to the full NLO result (full line) together with the 
separate contributions in Eq. (3.13) that are taken into account in moment space (dashed line) and 
in x space (dotted line), respectively. The latter is significant only for x — > 0. 



and 



dA 



(3.12) 



A=l-x 



It is straightforward to verify that, upon expanding Eq. (3.11) to 0(a s ), Eq. (3.8) repro- 
duces the full NLO result. The x — * 1 singular terms in — are taken into account 
through the Sudakov factor: they coincide with the 0(a s ) terms in Eq. (A. 4); the finite 
terms are incorporated into Mf^, and the remaining terms, which vanish at x = 1, appear 
explicitly in - ■ 

The explicit expressions for /ij,-(JV) appear in Appendix C. Having extracted the factor 
|.My| 2 which contains the constants at 0(a s ) in Eq. (3.8), ^{N) vanish at large N and 
the information needed to compute the 0{o? s IniV) term in Mjy T is fully contained in f(N) 
and the NNLL contributions to the exponent in Sud(m, N). 

Having matched the result at NLO only, non-singular terms at O(a^) are not under 
control. These terms depend on the details of the matching procedure. The particular 
choice we made gives preference to the moment-space treatment: as reflected in Fig. 7 the 
contribution of G2 is very small. 

To conclude this section we summarize the result. The Sudakov-resummed differential 



- 25 - 



rate, matched to the exact NLO result, takes the form: 
1 dT PT (x) r c+io ° 



dx 



dN /i/f pt , C F a s {m b ) 

2-KZ TT 



i,j = 1.. 

i < j 



->:i 



dA 



(3.13) 



A=l-a 



where the moments Mjy T are given by Eq. (3.11) and the normalization is fixed by r NLO , 
which is defined such that: 



r[6 — X sl ] Ei>Eo={l _ A)m/2 = r NLO x [0(A > 0) + R(A)} , 



(3.14) 



with 



1 NLO 



G ^ (m V t ;^i 2 (^w) 2 



Re{n} + 7 ^ ett ln 



(o) off , m 



where C^{jXb) = Cy ^ eff (//;,) + ^C^ eS (fif,). We emphasize that the contribution of the 
£ij terms, which is treated in x space, vanishes linearly with (1 — x) at large x and it is 
numerically very small in the experimentally relevant region of x. 

Finally, the corresponding formula for the total rate with an energy cut E 1 > E = 
(1 - A)m/2 is: 



1 v . f c+io ° dN (1 - A) 1 - N 

—T[b — » X sl ] E ^ >Eo = — { 



NLO 



2vri N —1 



M P N T + 



C F a s (m h ) 



TT 



^c^(A), (3.15) 



i,j = 1.. 



where c > 1. Note that at A — ► 1 (completely inclusive measurement) the inverse Mellin 
integral in Eq. (3.15) becomes trivial: it equals Mjv=i- 



4. Power corrections 

The perturbative analysis of decay spectra in the endpoint region [4, 16] exposes the pres- 
ence of three different scales: (1) hard, corresponding to virtual corrections with momenta 
0{m); (2) final-state dynamics of the jet with invariant mass squared 0(m 2 /N); and 
(3) soft radiation, with momenta 0(m/N), that determines the heavy-quark momentum 
distribution function. Naturally, there is also non-perturbative dynamics on each of these 
scales, which brings about power corrections. The largest power corrections are clearly 
those on the soft scale. Smaller power corrections, which nevertheless become increasingly 
important near the endpoint, are the ones corresponding to the jet-mass scale. Their com- 
bined effect is related to the appearance of resonances when the jet mass is sufficiently 
small. In the following we address power corrections on the soft scale only. 

4.1 Power corrections on the soft scale 

The classical approach to the quark distribution function [1,2] is based on the OPE. Ex- 
panding the non-local lightcone operator in Eq. (1.1) in powers of the lightcone separa- 
tion y~ one obtains local operators with a corresponding number of iD + derivatives. The 



- 26 - 



matrix elements of these operators in the meson define non-perturbative parameters f n , 
which, similarly to the non-local operator in Eq. (1.1), require renormalization. In the 
endpoint region all of these powers are relevant, and need to be summed. One arrives at 
the following result: 



(B(P)\ *(y)<f> y (p,yh+*(0) \B(P)) ~ exp 1 



2 ~ - ! M I 1 + E^V/M)' 

*> J L n=2 

(4.1) 

where the square brackets define the so called "shape function". Up to 0(1/N) corrections, 
the Mellin integral is the inverse of the Fourier integral in Eq. (1.1), so one gets the moment- 
space result by analytic continuation: iP + y~ — ► N. This leads to Eq. (1.4). 

A complementary point of view was presented in Ref. [16]. One begins by computing 
the quark distribution in an on-shell heavy quark, Eq. (1.2), which differs from the object 
of interest, Eq. (1.1), by the external states. Moments of this distribution are infrared safe: 
all logarithmic singularities cancel out between real and virtual corrections. However, the 
all-order result for these moments has infrared renormalons 11 , which indicate power-like 
sensitivity to large distance physics. As discussed in Sec. 2.2 above and in Refs. [16,30], 
the Sudakov factor Sn of Eq. (2.26) has infrared renormalons with ambiguities that scale 
as integer powers of NA/m: 



Sn = Sj/ x exp < 



J'=l 



where Kj represent 0(1) ambiguities. It was shown in detail in Ref. [16] that K\ has a 
special status: it is directly related to the leading renormalon ambiguity in the definition 
of the pole mass. This relation is realized in Eq. (1.4) by the fact that the normalization 
of the 0{N/M) term in the exponent is the difference between the meson mass and the 
quark pole mass. Higher-order terms in the exponential factor in Eq. (4.2) are reflected in 
Eq. (1.4) by T. 

Therefore, in principle the non-perturbative function T can be thought of in terms of 
the OPE or from the renormalon perspective. What is important in practice however, is the 
way the separation between the perturbative and non-perturbative regimes is implemented. 
Here we choose the Principal-Value prescription. We believe that this facilitates making 
maximal use of perturbation theory thus minimizing the role of power corrections. 

Based on Eq. (1.4), the physical photon-energy moments can be computed from their 
perturbative counterparts by 

M N = M P N T > PV x exp | - {N ~^ )Apv | x F((N-1)A/M); A PV = M-m PV , (4.3) 

with the matched expression for M^ T,PV given by Eq. (3.11), where the Sudakov factor of 
Eq. (2.17) is defined in the Principal- Value prescription. This implies of course the same 
prescription for A and for the higher-power ambiguities in T. In the numerical analysis that 



'See Ref. [30] for the large-/3o limit result. 
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follows we shall simply drop the unknown non-perturbative function T . Its parametrization 
would be worthwhile doing once stringent theoretical constrains are available or when 
experimental data become sufficiently precise. 

In this approximation we get the following DGE result for the differential spectrum in 
physical units: 



1 dT{E 1 ) _ M f c+io ° dN (2E^~ N 



dE^ 2 J c _ ioo 2iri \ M 



M N (4.4) 



Jc-ioo 2 ™ \m PV J 



PT,PV 

N ' 



where 12 Mjv in the first line is given by Eq. (4.3) with T = 1 whereas M^r T ' PV in the second 
is given by Eq. (3.11), where the Sudakov factor of Eq. (2.17) is defined in the Principal- 
Value prescription. While the two lines in Eq. (4.4) are trivially equivalent — the equality 
is violated by terms 0(A/M) — they reflect two different physical interpretations: in the 
first line the moments are the physical spectral moments which are free of any prescription 
dependence whereas in the second these are resummed perturbative moments in a given 
prescription. 

The equality of the two formulations in Eq. (4.4) indicates that the accuracy at which 
the meson mass is known is irrelevant for the calculation of the spectrum. On the other 
hand, no matter which of the two is chosen, an accurate determination of the pole mass in 
the Principal- Value prescription is crucial. 

4.2 Calculation of A PV 

In principle, A PV , or the Principal- Value pole mass m PV , can be determined from the 
relation with any other well-defined mass. A natural choice is the relation with the MS 
mass because it is reasonably well measured [48] 

mis(mMs) = 4.19 ± 0.05, GeV (4.5) 

and the corresponding perturbative expansion is known to the NNLO [45], while the nor- 
malization of the leading infrared renormalon is well under control — see Appendix B and 
Refs. [35-38]. 

Let us begin by defining m PV through the Principal Value of the standard Borel sum 
(see Appendix B): 

!^ = 1 + ^PV f dzB(z) e~ z / A , (4.6) 
mms Po Jo 

where A = /3oa s (mMs)/ 7r - Owing to the dominance of the first infrared renormalon, which 
is apparent already at low orders, and to the fact that the singularity structure of this 
renormalon is exactly known [34], it is possible to accurately estimate the normalization of 
the renormalon singularity [35-38] and thus construct a bilocal expansion of the form: 



b(z) = jr b n z n + 



n=0 



{l-2z) 1+ 2 5 



l+Yc k (l-2z) k 



k=i 



(4.7) 



12 For simplicity, we dropped the £ terms appearing in Eq. (3.13); while small, these terms are included 
in all the numerical results that follow. 
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where, as before, 5 = Pi/(3q. Here the first sum is known to the NNLO and we therefore 
use n max = 2. The coefficients in the second term depend only on the coefficients of the 
MS (3 function (see e.g. [35]) which sets fc max = 2. 

An accurate determination of the normalization constant q is essential for Eq. (4.7) 
to be useful for the calculation of m PV from Eq. (4.6). Clearly, also the convergence of 
the sums in Eq. (4.7) is relevant; but assuming that q is known, both series in Eq. (4.7) 
are free of any z = \ renormalon singularity 13 and thus their convergence is much better 
than that of the original m/m-gts expansion. Therefore, the accuracy at which m PV can be 
determined is comparable to the one at which ttims is known. 

According to Refs. [35-38] and Appendix B, q can be accurately computed. Using 
Eq. (B.20) we obtain: 

C F q(N f = 3)/vr = 0.560; C F q{N f = 4)/tt = 0.536. (4.8) 

This determination can be trusted within ~ 2 — 3%. Proceeding to evaluate m PV from 
Eq. (4.7) we obtain 



m f 



Til 

= 1.161 ±0.005; — ^ 

JV>=3 mMS 



= 1.164 ±0.005, (4.9) 

N S =A 



where we used the values in Eq. (4.8) for q, a s (m^ IS ) = 0.226, and estimated the error 
based on the NNLO contribution. Taking into account the effect of the finite charm mass 
we shall use m PV /mMg = 1.163 ± 0.005. Thus, based on Eq. (4.5) and the measured value 
of the mason mass, M = 5.279 GeV, we conclude that 

m PV = 4.874 ±0.050 GeV; A PV = M — m PV = 0.405 ± 0.050 GeV, (4.10) 

where the error is dominated by that of the short-distance mass determination, Eq. (4.5). 



5. Numerical results and comparison with data 

The state-of-the-art analysis of the total B — > X s j branching ratio is based on the NLO 
calculation [44,46,47]. According to Ref. [44] the integrated spectrum above raft/20 is: 

BR[ B _> X sl ] Ej>Eo=mb/20 ~ 3.73 • 10- 4 . (5.1) 

This prediction has ~ 10% accuracy. It is well known that a much higher cut on the 
photon energy is experimentally unavoidable. Realistic measurements in the B factories 
can be done with > Eq = 1.815 GeV [56], but higher cuts are advantageous and lead 
to smaller experimental errors. Moreover, it turns out to be useful to measure spectral 
moments over this limited energy domain. In this section we show that direct comparison 
of such measurements to theory is now possible. 

To describe the total rate we follow Ref. [44]. In this approach, instead of evaluating 
r NLO in Eq. (3.15) directly from its definition in Eq. (3.14), one uses the experimental value 

13 Higher renormalons are present in the first sum. In [36] the convergence of this sum was further 
improved using conformal mapping, which we do not use here. 
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for the B 



X c ev branching ratio avoiding explicit dependence on the b-quark pole mass. 



In addition, by an appropriate choice of mass scheme, Ref. [44] accounted for some charm- 
mass effects associated with charm quark loops in the O2 operator. These effects are 
formally NNLO in the strong coupling, but they are numerically important. Comparing 
the result of Ref. [44] with Eq = rrib/20, Eq. (5.1), to our calculation at the corresponding 
cut value A = 0.9, we fix the normalization of the branching ratio BRq: 



This formula will now be used to predict the dependence on A = 1 — 2£ , o/m PV . Throughout 
this analysis we shall use the central value 14 of BRo in Eq. (5.2). 

Recall that M PT ' in Eq. (5.2) is given by Eq. (3.11), where the DGE Sudakov fac- 
tor (2.17) is defined in the Principal- Value prescription. As discussed in Sec. 4.1 the use 
of the corresponding Principal- Value mass m PV in the definition of A guaranties consis- 
tency at the level of the "binding energy" ; it effectively accounts for the exponential factor 
in Eq. (4.3). On the other hand, additional power corrections o(^((N — 1)A/M) 2 ^ cor- 
responding to T in Eq. (4.3) are ignored. We are assuming that they are numerically 
small when using DGE with the Principal- Value prescription. Some evidence that this 
assumption is sensible is: 

• The DGE 'prediction is perturbatively stable: as shown in Fig. 8 the change from NLL- 
DGE to NNLL-DGE is moderate. This feature is best seen in moment space, see 
Fig. 4. Perturbative stability is of course only a necessary condition, not a sufficient 
one, as not all sources of power corrections can be probed by perturbative means. 

• The predicted spectrum smoothly extends beyond the perturbative endpoint and tends 
to zero for E y = (m + 0(A)) /2, close to the physical endpoint E 1 = M/2. In this 
respect the differences between the various models shown in Fig. 8 are less than 100 
MeV. This is quite remarkable given that in the close vicinity of the endpoint all power 
corrections on the soft scale are a priori important. On the other hand, this should not 
be taken as an indication that all power corrections are small: they may be large, yet 
conspire to cancel. Moreover, in the close vicinity of the endpoint the approximation 
based on keeping power corrections on the soft scale only breaks down: the non- 
perturbative structure of the jet becomes important; in our formulation resonances 
are simply ignored. 

These are, of course, qualitative arguments. In the long run data will hopefully allow 
quantifying these power corrections. 

The immediate advantage of the DGE calculation of the spectrum is the possibility 
to provide a reliable estimate of the partial branching fraction and higher moments with 

14 This is phenomenologically sensible because the sources of uncertainty in the distribution are different 
from those of the total branching fraction. 




i,j = 1..8 
i < 3 



with BR = 2.86 ± 0.29. 



(5.2) 
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Figure 8: The dependence of the differential (top) and integrated (bottom) B — > X s j spectrum 
as predicted by DGE on the model assumed for the Borel function of the soft Sudakov exponent 
(see Sec. 2.3). 



experimentally relevant cuts. To assess the remaining theoretical uncertainty in this pre- 
diction let us first compare the chosen model (c) for Bs(u) (Eq. (2.35)) to other models 
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Figure 9: The dependence of the differential (top) and integrated (bottom) B — > X s -f spectrum 
as predicted by DGE on the value of the strong coupling. 



for the soft Sudakov factor introduced in Sec. 2.3. Recall that the differences between 
these models are associated with the structure of the Borel function away from the origin, 
and are therefore equivalent in principle to the power corrections on the soft scale (m/N) 

15 Note that closer to the endpoint also the behavior of Bj(u) away from the origin becomes relevant. 
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E [GeV] 



Figure 10: The dependence of the differential (top) and integrated (bottom) B — > X s j spectrum 
as predicted by DGE on the assumed value for A. 



that would be parametrized by T. On the other hand, using the differences between the 
models as an error estimate is rather conservative 16 , since only model (c) is, by construc- 
tion, consistent with the large-order behavior of the Sudakov exponent. Fig. 8 shows the 

16 A more precise uncertainty estimate can be obtained by comparing different models that are consistent 
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differential and integrated spectra according to the various models for B$(u). Differences 
between models for the differential spectrum near the peak exceed 20%. For the integrated 
spectrum, however, this translates into less than 10% difference for practically any cut 
(excluding, of course, the immediate vicinity of the endpoint, Eq > 2.4 GeV). 



af s (M z ) 


< S (m PV ) 


Ajv /= 4 GeV 


A PV = M — m PY GeV 


0.113 


0.198 


0.265 


0.392 


0.1182 


0.216 


0.353 


0.405 


0.120 


0.222 


0.385 


0.455 



Table 1: Variation of the parameters entering the calculation of the spectrum with a s (Mz). 



Another theoretical-uncertainty estimate can be obtained by varying the input pa- 
rameters a s and to-ms within their error ranges. The sensitivity of the DGE spectrum to 
these parameters is higher than in typical short-distance calculations. The precise value 
of the coupling is important since a s is probed at low scales 0{m/N). The value of the 
short-distance mass mss directly influences the pole mass m PV (Sec. 4.2), which sets the 
scale for the spectrum in physical units, see Eq. (4.4). Table 1 details the values of the 
parameters entering the calculation for the central and the two extreme choices we made 
for af s (Mz). Fig. 9 shows the corresponding differential and integrated spectra. Fig. 10 
shows the spectrum for the central value of of the coupling, a s (m) = 0.216, while varying A 
by ±50 MeV. This corresponds, to a good approximation, to varying m^g within its error 
range. Note that while the variation in the coupling influences all the moments (the shape 
of the distribution gets modified) the variation of mas mainly affects the average energy: 
it generates a global shift of the distribution. 

A word of caution is due regarding the interpretation of the theoretical uncertainty 
reflected in the figures. As usual, uncertainty estimates are based on corrections that 
are known, e.g. degrading the NNLL-DGE to NLL-DGE. In this paper we have not 
parametrized non-perturbative power corrections, whose influence would increase towards 
the endpoint. Our analysis only reflects these effects through the different models for 
Bs(u). Such differences might not capture all possible power corrections. In this case, the 
uncertainty may be underestimated near the endpoint. 

Next, let us consider the average of the photon energy with a cut, namely 



7/E 7 >E 



f dE, 
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dTjE^ 
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dT(E 1 
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(5.3) 



with the NNLO expansion of Bs(u) as well as with Bs (1/2). Since the differences between the models 
introduced in Sec. 2.3 are of the same order of magnitude as other sources of uncertainty (e.g. the values 
of short-distance parameters) we postpone such detailed analysis to future work. 



^ 2.45 
> 

O 

uF 



at 

/ V 

UJ 

V 



2.35 



2.25 



Belle data 
NNLL-DGE (Model c) 

— A=0.355 GeV 
A=0.405 GeV 
- A=0.455 GeV 




2.2 2.3 
E [GeV] 



> 

at 
O 



0.06 



=0.05 



HI 

■ 

LU 



0.04 



0.03 



0.02 



0.01 



1.6 



_ I I I 


I I I I _ 

+ Belle data 




NNLL-DGE (Model c) - 




a s -0.198 




a s =0.216 - 




a s -0.222 






~ I I I I I I I I I I 


****** ' • ■ 

r i r i I i r i r I i i r i I i r i i I i i r i ~ 



1.7 



1.8 



1.9 



2.1 



2.2 
E 



2.3 
o [GeV] 



Figure 11: The first two truncated moments of the B — ► X s j spectrum: (E~ ( ) (top) and 
(((Ej) - E~ ( ) 2 ) (bottom), as a function of the minimum photon energy cut Eq, as calculated by 
DGE using model (c) of Eq. (2.35) while varying a s and m^s within their error ranges. In each 
plot we show a few curves representing the dominant source of uncertainty considering the various 
models for Bs(u) and the values of a s and Toms- The DGE result is compared to the available data 
from Belle [56] with a cut E = 1.815 GeV. 



and, similarly, higher truncated moments: 
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Figure 12: DGE results for the first ten truncated moments of the B — ► X s ~f spectrum, 
(((-Ey) —E y ) n ), with a cut E 1 > Eq = 1.8f5GeV. All moments are presented in energy units 
(GeV) by taking the corresponding n-th root. The n — 2 result is compared with Belle data, where 
the internal and external shaded regions represent statistical and statistical plus systematic errors, 
respectively. 



Note that these moments differ from the standard Mellin moments discussed so far in two 
respects: first, they are defined over a limited photon-energy range, and second, these are 
"central moments" in the sense that they depend of the difference between the photon 
energy and the average energy in this range. 

For sufficiently high Eq cuts and low enough n these moments are accessible experi- 
mentally. Results for the average energy, Eq. (5.3), and the variance (n = 2 in Eq. (5.4)) 
were recently published by Belle [56,57] with Eq = 1.815 GeV. In Fig. 11 we compare this 
experimental result to our calculation while providing additional theoretical predictions for 
the dependence on the cut value. We find very good agreement. The experimental error is 
somewhat larger than the theoretical one. Note that the dominant source of uncertainty 
in the theoretical prediction for the average energy is the value of ttt-msj affecting the cal- 
culation of the pole mass m PV or A PV . In the future it may even be possible to get an 
accurate measurement of rass from experimental data for {E^) E >Eq . On the other hand, 
the dominant source of uncertainty in the variance is the value of the strong coupling. 
More interesting comparison can hopefully be done in the future with experimental data 
having higher cuts, where the systematic experimental errors are expected to be smaller. 
The possibility to make comparison between theory and data as a function of the cut has 
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an important added value: clearly, at sufficiently high Eq values the theoretical prediction, 
which is lacking any non-perturbative corrections, will fail. 

Finally, in Fig. 12 we show the theoretical prediction for the first ten moments. Going 
to higher moments the finer details of the shape become important. The expectation is, 
of course, that power corrections will become increasingly important at high n. Having 
excluded power corrections, the dominant source of uncertainty in Fig. 12 is the value of a s . 
Since for higher moments the coupling is probed at lower scales ~ 0(m/n), this uncertainty 
gradually increases with n. Following Belle [56] we have chosen E-y > Eq = 1.815 GeV. 
Fig. 12 shows the comparison with Belle data for the variance; there is no data to compare 
with for n > 3. 

6. Conclusions 

Inclusive -B-decay spectra present a special challenge to theory: because of their inclusive 
nature the details of the hadronic wave function are largely irrelevant, but on the other 
hand, the soft scales involved prohibit a straightforward perturbative approach. 

• The clearest manifestation of the non-perturbative nature of the B — ► X s j spectrum 
is that any fixed-order perturbative result has support only for E 1 < m/2 while the 
physical spectrum extends to E^ = M/2. The gap between the quark pole mass and 
the meson mass is related to non-perturbative dynamic of the bound state. 

• Moreover, the perturbative endpoint region is characterized by parametrically large 
corrections. The conventional approach to Sudakov resummation, which relies on the 
dominance of logarithmic corrections in the large- A" limit where a s {m) In A is fixed, 
fails because of large running-coupling effects. The breakdown of the perturbative 
expansion with increasing logarithmic accuracy is demonstrated in Fig. 2 (see also 
Fig. 2 in [16]). 

In this paper we showed that DGE provides a solution to both these problems. Al- 
though it is based on resummation of perturbation theory the DGE-resummed spec- 
trum does not have perturbative support. By taking the Principal- Value prescription 
in moment space this inherent limitation of any fixed-order approximation is removed. 
The spectrum smoothly extends to the non-perturbative regime and tends to zero for 
Ej = (m + 0(A)) /2, qualitatively as expected in the meson decay (see Fig. 6). 

Of course, the DGE spectrum does not have correct non-perturbative support proper- 
ties: it does not strictly vanish for Ey > M/2. Moreover, is not necessarily positive definite 
for any E 7 < M/2 (see e.g. Fig. 8). Indeed, the close vicinity of the endpoint is beyond 
theoretical control as it involves an infinite number of non-perturbative power corrections 
on the soft scale m/N, and eventually also the non-perturbative structure of the jet. 

In this work we presented perturbative predictions and refrained from making any 
parametrization of non-perturbative power corrections. It is conceivable that there exist a 
region where a few non-perturbative corrections on the soft scale are phenomenologically 
relevant, while corrections on the jet-mass scale can be ignored. Once experimental data 
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are precise enough, power corrections on the soft scale will be worthwhile parametrizing be- 
cause of their universal nature: they are related exclusively to the quark distribution in the 
meson. Therefore, fixing their magnitude can be instrumental to precision measurements 
of V u b from charmless semileptonic B decays. The DGE spectrum presented here provides 
a baseline for systematic parametrization of such non-perturbative power corrections. Per- 
turbative stability is imperative in this respect. The remarkable stability of the DGE result 
(compare Fig. 4 to Fig. 1) is largely thanks to the resummation of running-coupling effects. 

It should be emphasized that the difference between DGE and any fixed-logarithmic- 
accuracy approach is a qualitative one. Recall, that an accurate value of m PV (or A PV ) is 
essential for the calculation of the photon-energy spectrum in physical units (see Sec. 4). As 
opposed to the general concept of the pole mass, m PV has no linear 0(A) ambiguity. It can 
be determined with roughly the same accuracy as short-distance masses (~ 50 MeV). It is 
only because the Sudakov exponent was defined using the Principal- Value prescription that 
m PV becomes the relevant mass. This is a manifestation of the cancellation of renormalon 
ambiguities between the quark-distribution subprocess in the Sudakov factor and the pole 
mass [16]. In a fixed-logarithmic-accuracy approach the renormalon in the Sudakov factor 
is hidden — it is realized through the divergence of the expansion in Eq. (2.11) — so the 
choice of mass scheme when computing the spectrum in physical units becomes ad hoc. 

Experimentally-favorable observables are moments of the photon-energy spectrum 
defined over a limited range, E 1 > Eq, Eqs. (5.3) and (5.4). Having obtained a stable 
prediction in Mellin space over a wide domain in the complex N plane, we can reliably 
compute the truncated moments over a range of experimentally-relevant cut values. We 
showed that our calculation nicely agrees with the recent data from Belle [56]. Furthermore, 
a great variety of possibilities for comparison between data and theory is now open, e.g. 
studying the dependence on the cut value and higher truncated moments. In this way 
more detailed information on the distribution can be systematically extracted. This is, in 
particular, imperative to quantifying power corrections. 

The qualitative success of the DGE calculation presented here, which does not involve 
any non-perturbative parameters, and the good agreement with the available data from 
Belle, indicate that additional power corrections in this framework are indeed small. This 
approach is therefore promising for extracting V u b from charmless semileptonic B decays. 

Note added 

A few weeks after the submission of this paper, the BaBar collaboration has presented [58] 
new preliminary data for the partial branching ratio, as well as the average energy and the 
variance as a function of the cut. Comparison of the predictions in Figs. 10 and 11 above 
with these measurements can be found in Ref. [59]. 

Acknowledgements 

EG wishes to thank Vladimir Braun, Gregory Korchemsky and Arkady Vainshtein for very 
interesting discussions and the KITP at UC Santa Barbara for hospitality while some of 



- 38 - 



this work was done. JRA acknowledges the support of PPARC (postdoctoral fellowship 
PPA/P/S/2003/00281). The work of EG is supported by a Marie Curie individual fellow- 
ship, contract number HPMF-CT-2002-02112. 



A. Singular terms at NNLO 

Let us expand the Sudakov exponent in Eq. (2.11) to O (a 2 ). We get: 



Sud(m, N) = exp 
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To get M^ T '° 7 in Eq. (2.1) one needs the matching coefficient, which is currently 
known only to 0(a s ). It is given by [16]: 
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Let us note that having fixed the anomalous dimensions as above, the only missing ingredi- 
ent for Sudakov resummation of decay spectra with NNLL accuracy is the ^-independent 
term at 0(a 2 s ). 

For easy comparison with future two-loop calculations of the decay process, we sum- 
marize below the log-enhanced terms to order a 2 . Expanding the exponent and taking into 
account the constant terms at order a s in Eq. (A.2) we obtain the following log-enhanced 
terms at O {a 2 ): 
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Note that the normalization of the rate influences not only constant terms at 0(a s ) but 
also the contributions to the In 2 N and the In 1 TV terms at 0(a 2 s ) that are proportional 
to Cp. Here we present the expansion of Mj^' 7 , where the N = 1 moment is exactly 
1 by definition. It is straightforward to convert this expansion to other normalization 
conventions (Eqs. (C.2) and (3.14) may be useful). 

Converting Eq. (A. 3) to x space, the x — > 1 singular terms in the b — ► X s ~f distribution 
(from Ot) up to 0{a 2 s ) are: 
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We checked that the terms that are leading in [3q in Eq. (A.4) agree with the large-x 
expansion of Ref. [8] . 

In Eq. (A.4) { } + stand for the standard "+" prescription, i.e. 



{p(z)} + = p(z) - 5(1 - z) / p(z)dz, 

Jo 
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so, writing 
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the convolution with a smooth test function f(x) (representing non-perturbative correc- 
tions) takes the form: 
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where both terms are well defined. In moment space, this is equivalent to multiplying 
the moments of the perturbative distribution Mjy T by the moments on the test function: 
F N = f 1 dzz N - 1 f(z). 
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B. The u — \ renormalon in the pole mass 



In this Appendix we determine the normalization of the u= \ renormalon in the pole mass. 
This result is used for the calculation of the pole mass in the Principal- Value prescription 
in Sec. 4.2 and for the comparison with the soft Sudakov exponent in Sec. 2.3. 

Remarkably, the structure of the leading renormalon ambiguity in the pole mass is 
known exactly. Ref. [34] has shown that, owing to the vanishing of the anomalous dimen- 
sion of the h v h v operator, the linear ultraviolet divergence in the self energy of a static 
quark has a very simple structure: it is just proportional to the QCD scale A without 
any logarithmic corrections. Consequently, the imaginary part associated with the u = \ 
infrared renormalon in the pole mass has the same property: 

Im{m} = const x A. (B.l) 

This means that, up to an overall normalization constant, the large-order behavior of the 
relation between the pole mass and any short-distance mass (which is induced by this 
renormalon) depends only on the coefficients of the f3 function. 

Owing to the dominance of the u = \ renormalon, which sets in already at low orders, 
the normalization constant can be computed from the perturbative expansion with accuracy 
of a few percent [35-38] . 

The renormalon in the standard Borel representation 

Let us first briefly review the standard analysis (more details can be found in [49]). To 
make use of Eq. (B.l) one first integrates the renormalization-group equation, 
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where a(p) = /3oa s (/j,)/ir and 5u = Pk/Po +1 , writing A as: 



A // oxp \ ^- - - hi(ad) 



k=l 



(B.3) 



ci 

C3 



where are specific combinations of the coefficients of the (5 function; for example: 
- 6 2 + 6 2 

2 63 + 4 6 2 6 - 2 5 s + 5 2 2 - 2 5 2 5 2 + 5 4 
- 864 + 66263 - 5 2 3 + 85 2 2 + (16S 3 - 126 2 2 )6 



1 r 
2 

1 " 

8 - 
1 

48 



+(3£ 2 2 - 66 3 - 24 6 2 ) 6 2 + 18 6 3 6 2 + (8 - 3<5 2 ) <5 4 - 6<5 5 



(B.4) 



Then, writing the perturbative relation between the pole mass and the MS mass in the 
standard Borel representation (where the Borel variable z is conjugate to a), 
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the singular part of B(z) near z = \ can be explicitly written using Eqs. (B.l) and (B.3): 
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To see this one inserts Eq. (B.6) into Eq. (B.5) and takes the imaginary part of the integral, 
getting: 
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In this calculation one uses the following formula: 
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which is derived assuming that zq has an infinitesimally small positive imaginary part. 
Eq. (B.8) clearly has exactly the same dependence on the coupling as Eq. (B.3), so the 
requirement of Eq. (B.l) is satisfied. 

The pole mass (and its imaginary part) is renormalization-scheme invariant. Several 
authors have computed the normalization constant q using the perturbative expansion of 
"V^ms^ms) i n the MS scheme [35-38]. This perturbative relation is available to 0(a^) 
— see Eq. (10) in [41] (the MS mass anomalous dimension has been computed before [42]). 
Knowing that z = \ is the nearest singularity to the origin and using Eq. (B.6) one can 
extract the residue by first multiplying [39, 40] the available n-th order approximation 
of B(z) by (1 — 2z) 1+ 2 <5 , then expanding in powers of z and truncating at order n and 
finally substituting z = \. We performed this calculation and we agree with previous 
results. The residue is shown in Fig. 13 as a function of (5q. For example, we obtained: 
Cpq(Nf — ► oo) /it ~ 0.935 and Cpq(Nf = 4)/7r ~ 0.551. We estimate the error on this 
determination as ~ 3%. This estimate is based on comparison with a similar calculation 
using the scheme invariant Borel transform (see below) and on the one point where the 
exact value of the normalization constant has been computed [27], the large-/3o limit. There 
B(z) is 



B{z) 



largeft) 



3 5 2 (1 - z)V(z)Y(l - 2z) 
regular + . 
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where the regular terms are related with the renormalization of the MS mass. Thus the 
exact value is Cpq{Nf — ► oo)/7r = C^ee/V ~ 0.9766. 

To clarify the relation between the normalization constant in the scheme-invariant 
Borel transform and in the standard one, let us first consider the latter in different renor- 
malization schemes. First, let us stress that having fixed the definition of A in Eq. (B.3), 
the scheme (and thus the value of the coupling at the scale fi) depend on all the /^-function 
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Figure 13: The normalization constant of the u 



i renormalon of the pole mass (Cf^/V) as a 



function of 1//3q- The exact result of the large-/?o limit is denoted by a cross. Calculations relying 
on the perturbative expansion of m/m-^(m-^): standard Borel transform of the mass ratio in MS 
(dotdashes); scheme-invariant Borel transform of the mass ratio (full line); the latter optimized 
using Eq. (B.20) (plus signs). Squares represent the calculation relying on the scheme-invariant 
Borel transform of the soft Sudakov exponent and on the cancellation of the ambiguity (i.e. on 
Eq. (2.30)). This calculation uses an optimization procedure similar to the one of Eq. (B.20). 



coefficients (5k with k > 2. We note that while the sum in Eq. (B.6) converges, the signi- 
ficance of subleading terms there strongly depends on the scheme. In particular, we find 
that in a scheme (RS) where the /3-function coefficients are characterized by geometrical 
progression, da RS /dln/x 2 = — o| s / (1 — 5a RS ), all the coefficients vanish. On the other 
hand in the 't Hooft scheme of Eq. (2.18), one obtains cjj Hooft = (S/2) k /kl, so the singular 
part of the Borel function sums up into: 
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so 17 the imaginary part in Eq. (B.8) sums up into: 

where the 't Hooft coupling A is evaluated at mis(™is)' The same result can be obtained 
by substituting Eq. (B.ll) into Eq. (B.5) (where a is replaced by the 't Hooft coupling A) 
and evaluating the integral using the relation: 

where T(u) and the relation between A and In ^ 2 /A 2 are given in Eq. (2.18). The imaginary 
part of the l.h.s. is simply — ir times the residue at u = \. This exercise also shows that the 
scheme invariant formulation is most natural for the problem at hand: the singular part of 
the corresponding Borel function is just a simple pole. 

The renormalon in the scheme invariant Borel representation 

Let us turn now to consider the pole-mass renormalon in the scheme-invariant formulation 
of the Borel transform. The mass ratio is expressed as 

C [°° ( A 2 \ u 
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where as in Eq. (2.18), T{u) corresponds to the 't Hooft scheme. 
Taking the imaginary part of Eq. (B.14) we obtain: 
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where e > is infinitesimally small; the two expressions were obtained based on the residue 
theorem and Eq. (B.13), respectively. In the second expression A(jx) is defined by 

2 = 2 + 4 (B.18) 



7 Note that using these coefficients in Eq. (B.3) one recovers the definition of A in Eq. (2.18). 
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or, equivalently, In (/U 2 /A 2 ) = 1/A - 5 + <51n(<L4). Comparing Eq. (B.15) with Eq. (B.12) 
we find that 

q = -2Re~ 5/2 . (B.19) 

The approximate calculation of the residue R based on the perturbative expansion 
proceeds similarly to the standard Borel transform: one multiplies by — u) expands 
and uses the truncated series to compute the value at u = \. The result is shown in as 
a function of 1/(3q in Fig. 13, where the comparison with the standard Borel transform is 
available. One can repeat the calculation by first multiplying the truncated series by an 
arbitrary function f(u) (which has a Taylor expansion at u = that converges at least for 
\u\ < tj), then taking the limit, and finally dividing by the exact value f{\), 

R opt = -L ton _ B(u)f(u). (B.20) 

We have chosen the set of functions: f(u) = exp(rn) where r is arbitrary. A reliable 
approximation should be independent of variations of the function f(u). Therefore, the 
procedure can be improved by finding a saddle point with respect to r. The result of this 
optimized calculation of R are shown in Fig. 13 as well. It turns out that the result is quite 
stable so the difference between the optimized calculation and the simple one is small, for 
example, for Cpq(Nf = 4)/-7r we get 0.5355 and 0.5363 in the two, respectively. 

In conclusion, the normalization of the u = \ renormalon in the pole mass can be 
determined from the perturbative relation with the MS mass within a few percent. The 
difference between the determination using the scheme-invariant Borel representation and 
the one using the standard Borel representation is somewhat larger than the variation of 
the results within each of the two approaches. 



C. NLO results for b — > s^g and the matching procedure 

We begin with the well-known expressions for the b — ► s^g contributions [45]: 
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where 
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and where, following [44], we use z = m^ s (n)/rrib ^ 0.22 ± 0.04. 

As explained in Sec. 3 our matching procedure involves splitting some of these func- 
tions, 4>ij{A) = i]ij(A)+£ t ij(A), such that rjij(A) dominates the small-A limit. We require: 
^•(A) = 0(A 2 ). 

For 077, which is conveniently transformed to moment space we simply define 7777(A) = 
77 (A) and £ 77 (A) = 0. It then follows from Eq. (3.12) that 



Htt(N) = 



(*(7V) + 7E ) 



1 1 

1 + N + N 



1 1 



9 

+ — + 



+ 



2N + 1 N + 2 2N (AT + 1) 2 N 2 



Note that in Eq. (3.11) the full O7 contribution is reproduced by f(N) + /j,ji(N) plus the 
single and double logs from the expansion of Sud(m,A) (these are the 0{a s ) terms in 
Eq. (A.3)). Note also that f(N) + ^77 (AT) differs from C?7 of Eq. (A.2) by the constant 
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As explained in the text, the treatment of the other contributions involves splitting 
the real-emission terms into two. For all the contributions associated with the operators 
0\ and O2 and their interference with O7 and 8 (this includes 022, 027, 0n, 012 017, 018 
and 028) we define 



"*3 



dA 



A=0 



and write: 



Vij(A) = 4>ij A, 

iij{A) = 4>ij{A) — 4>ij A. 

Since the derivative of %(A) is a constant, the moments in Eq. (3.12) are pure 1/N terms: 
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Finally, the contribution of 088 it taken into account by extracting the leading term 
at A — ► (which makes the corresponding contribution to the differential rate singular) 
and incorporating it in moment space: 
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As follows from their definition all Hij(N) vanish at large N. 

Finally, let us summarize the result for the first moment which enters the expression 
for the total rate (Eq. (3.15) with A — ► 1): 
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The remaining /iij(N = 1) values can be obtained from the above using the relations in 
the last line in Eq. (C.l). 
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